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ff^ ■ We suggest two methods (based on flavor SU{3) symmetry) to 

Q ■ determine the CKM angle 7 using the decays B^, B^ — > Ktj (r/') and 

Qn ' Bg ^ TTT] [r]'), respectively. Rescattering efl'ects are partly included 

- we neglect annihilation amplitudes, but do not assume any other 
p 1' relation between the SU{3) invariant amplitudes. We use the fact 

p ' that the amplitude (including the Electroweak Penguin contribution) 

for Bii, B^ — > ttK with final state I (isospin) = 3/2 is known as a 
function of 7 from the decay rate B'^ -^ tt^tt'^ . 
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1 Introduction 

The principal aim of the B physics experimental programs is to measure the 
angles (denoted by a, (3 and 7) of the triangle representing the unitarity 
relation: V^^Xtd + V*^Vcd + KT^Kd = 0, where V is the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. The idea is to overdetermine the angles of this 
triangle and thus test the CKM paradigm of CP violation. 

Methods to determine 7 (= Arg (— K*6Kd/Kfe^cd)) including the Elec- 
troweak Penguin (EWP) diagram contribution have been suggested [^ |], ^ 
§, ^. To (over) determine 7 to test the CKM theory of CP violation, we mea- 
sure 7 using different techniques. It is thus useful to have new methods to 
determine 7. With this motivation, in this letter, we give two new methods 
to determine 7 using time integrated rates for the decays B^,, B^ -^ Kr] [tj') 
(Method 1) and B^ -^ ttt] [t]') (Method 2). As in all the other methods 
mentioned above, flavor SU{3) symmetry is used in both the methods. 

We will write the decay amplitudes in terms of flavor SU{3) invariant 
amplitudes 0. These are denoted by C3 ' , Cg ' , C^g , ^3' and A^i, 
and correspond to the 5 linearly independent ways of forming flavor SU{3) 
singlets from the initial meson Bi, the two final state mesons belonging to the 
flavor SU{3) octet and the effective weak Hamiltonian which transforms as 
a 3 X 3 X 3. T and P denote the parts of these amplitudes generated by tree 
level and penguin operators respectively. These invariant amplitudes include 
soft final state rescattering effects. Some of the methods to determine 7 [^], 
including the ones which use the decays B^ -^ Kr] {rj') [^ |[ and the decay 
Bs — > 7r?7 {t]') |jl| neglect rescattering effects. In particular, these methods 
assume that the decay amplitude B^ -^ n^K^ has no weak phase e*''' from the 
tree level operators. In the language of the SU{3) invariant amplitudes, this 
is equivalent to assuming that the annihilation amphtudes | Ai are suppressed 
by /b/'^b and a combination of the SU{3) invariant amplitudes, Cj — Cj — 



^Annihilation amplitudes are the ones in which the index i of 5, is contracted directly 
with the Hamiltonian. 



C^ is zero (both of which are vahd in the absence of significant rescattering 
effects). Rescattering effects can enhance the annihilation contributions and 
lead to significant Cj — Cj — Cf^ . In this letter, we neglect annihilation 
contributions but do not assume any relation between Cj, Cj and C^ or 
the other S'f/(3) invariant amplitudes. Thus, rescattering effcts are partly 
included. 

The decay amplitudes for B^ — > "kK can be written as H] 



-A{B,^.-K^) = Xi'\Cj + Cj + 3C^,)+Y.KHc[, + C^,, + ^Ci 
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q 
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Here, A^'^'-' = V^^Vg^/ (g = u,c,t and q' = d, s) and C^, A^ denote the penguin 
amplitudes due to quark q running in the loop. Using the unitarity of the 
CKM matrix, i.e., A^^ = -Xi'^ - Xi'\ we get 

Xl'^Cr + E K'K, = >^CJ - A(^)Cf , (2) 

where Cf = Cj — Cf^ + Cf^ and Cf = Cf^ — Cf^. A similar notation is used 
for Aj and A[. Henceforth, we will write the decay amplitudes using this 
notation. 

The decay amplitudes for B^ — > ttK are [^ 

AiB^^n^K') = Xl^^Cj-Xi^^C[-Xi^^Cj + X^:^C^ 

~X^ Ci^ + A^ Ci5 + 3A„ y4]^5 — 3A^ A^^,, 
-V2A{B^^n'K^) = AWC3^ - A(^)C3^ - A^CJ + A^Cr 

+7X(:^Cl - 7X^C^ + 3Ai^)if, - 3AWAf,. 

(3) 



Using Eqns.(P and (|^), we get an expression for A3/2, the amplitude for B^ 
Bd -^ ttK with final state / (isospin) = 3/2, 
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n'K- 



Oj^O) 



A{B+ ^n+K^) + V2AiB^ 
= A{Bd^n-K+) + V2AiBd 
= -8 (AWcf, - A^Cf,) . 

The decay amplitude for B^ -^ ti^tt^ is 

- V2AiB^ - TT+vr") = 8 (Af C?, - \^^C[,). 



(4) 



(5) 



The QCD penguin diagram (which is A/ = 1/2) does not contribute to this 
decay since this decay has the transition A J = 3/2. So, Cf^ is the EWP 
contribution. Neubert, Rosner [Bi showed that 






tvic: — l^n 



(6) 



where Kg = (cg^^ + Cio,g)/(ci + C2) is the ratio of Wilson coefficients (WC's) 
of the EWP operators (with quark q running in the loop) and the tree level 
operators in the effective Hamiltonian so that 
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Af [l + l-u- l^t 
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The top quark EWP diagram with Z exchange is enhanced by iri^/m\ and 
so Kt > Ku^c giving 



^/2A{B^ 



n TT 



0\ 



clxi'^ 



15^*u 



2 



AW 3 

--K 



y(d)2 



where k = Kt. Since 3/2 k ~ 2% and jA^'^^l ~ lA^"^^], we get 



V2A{B' 



tt+ttO) 



cl\>^f'V 
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(8) 



(9) 



in the Wolfenstein parametrization and setting the strong phase of Ci^ to 
zero, i.e., the EWP contribution is ~ 0(2%) and can thus be neglected in 



the decay amplitude 5+ -^ 7r°7r+. Q Thus, C^ can be determined directly 
from the decay rate B'^ -^ 7r+7r°. 

Similarly, the expression for A3/2 (Eqn.(^) simplifies to 



A3/2 ^ -8 Cl (\Xl:Y^ - mL^ (10) 



A3/2~-«0'{5(^|Ar|e"-|Ar' 
so that 



1^3/2! = 8 Cf5|Ai'^V(l + 61^ - 26ewcos^), (11) 

where 6ew is given by |A^''^|/|A^''-'| 3/2 k ~ 0(1), i.e., in this case, due to the 
CKM factors, the EWP contribution oc A^*-* «s important. Thus, knowing C^ 
from the decay rate B^ -^ tt+tt^ and 7, we can determine A^/2 and conversely 
7 can be determined if the (magnitude) of A^/2 is known using some other 
method and the decay rate B^ -^ tt^tt'^ is measured. In using this relation, 
it is crucial that the parameter 6ew is calculable. 

In the analysis up to now, annihilation contributions are included. 

2 Method 1 

The decay amplitudes for Bd, B^ decays to rjs = 1/VQ (2 ss — uu — dd\ and 
?7i = l/-\/3 {ss + uu + dd\ can be written as 

-9 {XI'^CI - A(^)Cf,) + 3 (Ai^)if, - A^Af,) , 
V3A (5+ -> K^vi) = 2 {xl^^Ci - Xi^^C[) + Al^)Cj - Xi^^C^ 

+3 (AWcf, - X[^^C[,) + 6 {xi^^Al - A(^)ylf,) 
+3 {XI^^EJ - X[^^E[) + 3 {xi^W^ - Xi^^i) 
+9 {XI:^DI - A(^)Df,) , (12) 



^This EWP contribution can actually be included g], but we neglect it for simplicity 
here. 



^15 y ; 



—5 (^A„ C]^5 — A^ C15J — (^A^ A]^5 — A^ A^ 
VSA {B, -. K\) = 2 {Xi^^cj - Xi^^Ci) - Xl^^Cj + X^f^C^ 

- {^i'^Cl - A(^)Cf,) - 2 (A^if, - Xi^^A{ 

+3 (A(^)EJ - AWE3O - 3 {^L'^DI - ^i'^D^ 
-3 (Ai^)D?; - A^Df,) , (13) 

where E3, D^ and D15 are the amphtudes which contribute only to B meson 
decays to a final state involving r/i [0. ^ Dq and Di^ are annihilation 
amplitudes. We assume that the mass eigenstates i] and tj' are given by the 
the canonical mixing []: 

2V2 1 

V = -^Vs--^Vi 

= —= (ss — uu — dd\ , 



1 , 2^2 

ri = -^vs + ^m 

= -^{2 SS + UU + d(tj . (14) 

This mixing is consistent with the present data 0. Then, the decay ampli- 
tudes for B^ , Bd decays to 77 and rj' are 

-V?>A{b^^ K^r^) = X^^Cl - X^C^ + 7 (A^C?, - X\^^C^,) 

+3 {XI'^DI - Xi^^D^,) , 
V6A{b^^ K^i) = 3 {XI:^CJ - Xi^^C[) + X(:^CJ - X^^C^ 



^The notation D'^ , Dp is similar to C*^, Cp (Eqn. 
^Both the methods can be easily modified in the case of a general mixing angle, provided 
the mixing angle is known. 



+12 (Ai^)Df, - A^Df,) , (15) 



'15 



- v^ ^ (5. - ir%) = -A(^)C,^ + Xi^^C[ + 3 (a(^)C?; - A^Cf, 

+ (Ai^)^3^ - X[^^E[) - (Ai^)Dj - A(^)D, 

- {>^l:'DL - Ai^^/^fs) , 

V6A{b,^ KW) = 3 (a(^)C3^ - Xi^^C,^) - X(:^CJ + A(^)Cf 

-3 {Xi^^cl - Xi^^C[,) - 3 (A^if, - A(^)A 
+4 {XI:^EJ - Xi^^E^) - 4 {xi^W^ - A^D, 
-4 (XI:^DL - A^Df,) . (16) 
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From Eqns.(D, ([T|) and (|l|), we get the relation g 

Ve^ (5+ -^ K+ris) = ^VSA (5+ -^ K+ii) + ^Ve^ (5+ -^ K+?]' 

= 2A{B+^ 71+ K^) + V2A (B+ -^ n°K- 



(17) 

As mentioned earlier, the magnitude of A3/2 is known as a function of 7 in 
terms of the B+ -^ 7r°7r+ decay rate (Eqn.(pTD). This, for a given 7, we can 
construct the two triangles formed by B^j -^ ttK, A3/2 and B+ -^ nK, A3/2 
(corresponding to Eqns.(^ and A's DEB and ADB respectively in Fig]^) 
and the quadrangle formed by B+ -^ 71K, K+rj, K+rj' corresponding to 
Eqn.(|l^) {AD EC of Fig.[l|). Thus, we know the phases (in the convention 
where the phase of C^ is zero) of the decay amplitudes B+ -^ ttK, riK+, 

®There are two discrete ambiguities in the construction of Fig.H: A's ADB and DEB 
can be on the same side of the common base DB and similarly in the quadrangle ADFC 
the vertices A and F can be on the same side of the diagonal DC. 
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1: The polygon construction for Method 1. 
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1^6 \A{B+ ^ K+i)l GC = 3 FG = V3 \A{B+ ^ K+r])l DE 

\A{Bd -^ n-K+) \. Given 7 and the B+ 



0^0 



BE 



V2 \A{Bd^7r 

tt+tt'^ decay rate we know DB 

VS\A{Ba^K^r])\ is EG. 



A^/2\ from Eqn.fpTI). The prediction for 



K^Tj' and Bd -^ tcK as a function of 7 form this construction (the magnitudes 
are, of course, known from the measurement of the decay rates). 

We have included the annihilation contributions up to now. 

If we neglect the annihilation amplitudes, Ai^, Dq and D15, using Eqns.(|I]), 



(|T5D and ([T6|), we get the relations 



VqA (Bd -^ K\ 



-V2A[Bd^ix'K') + -VQ{B+ 


-. K^v') 


+ \V^A[B+^K+r^), 


(18) 


- V2A [Bd ^ 7r°ir°) + ^Vq (i?+ - 


. K^V') 


+^VzA[b+ ^ K+^i) . 


(19) 



Thus, we can predict the decay amplitudes Bd — ^ K^r], K^rj' as a function of 
7 since, as mentioned above, we know the magnitudes and phases (the latter 
from Fig.|lD of all the other decay amphtudes in Eqns.(|18D and ([T9|). In fact, 
the decay amplitude Bd — > K^r] is shown as EG in Fig.0. Once the decay 
rate Bd —>■ K^rj or K^rj' is measured, 7 can be determined. Thus, 7 can be 
determined (up to a four-fold discrete ambiguity) by the measurement of the 
decay rates for 8 modes - B^ —>■ tt^7t~^, ttK, K^rj, K~^ri', Bd — > tcK, K^t] (or 

3 Method 2 

This method is based on the method of Gronau et al. [jl|. In |l|, the an- 
nihilation amplitudes are neglected and the relation Cj — Cj — C^ = is 
assumed. In other words, rescattering effects are neglected so that the am- 
plitude for the decay B^ -^ K^t^^ has no weak phase e^"' . We neglect the 
annihilation amplitudes, but do not assume any relation between the C's. 
Thus, we include partly the rescattering effects. The decay amplitudes for 



Bs -^ 7r°?78, ?7i are 

+4 (Al^^ifs - Ai^^Afs) , (20) 



V6A{B, - vr^O = 2 (A(^)CJ - A(^)Q^) - 4 {xi^^Cf, - Xi^^C[,) 

-8 (Ai^^if^ - Ai^Mf5) + 6AWdJ - GAW^r 
-12Ai^)Df5 + m^Df^. (21) 

With the canonical mixing (Eqn.|l^), we get 

V6A{B, - vr%) = 2 (A(^)CJ - A^Ce^) - 4 (A^C?; - A(^)Cf,) 

+8 (Ai^)ir5 - Ai^^f^) - 2 (A^DJ - A(^)Df ) 
+4 (A(^)Df, - A(^)Df,) , (22) 



^(5,^ A') = v^(5. - A) - 2v^ [2 (AWifs - A^^^Af^^ 

- Al^)^6^ + A^Z^e^ + 2 (A(^)I)?; - A(^)Df,)l . (23) 



Neglecting the annihilation ampltudes, Ai^, Dq and D15, from Eqns.([l|), (|^), 
and ( [23| ) we get the relations |I[]: 



VQA (b, -^ 7r%) ^ -A (b+ -^ K\+) + V2A (Bd -^ K' 



o„o 



= V2A [B+ -^ i^+7r°) -A(Bd^ K+n-) , (24) 
Ve^ (b, -^ 7r%) ^ V3^ (b, -^ 7T%') . (25) 



From Eqns. ([2^) and ( P^D we see that the decay amplitudes B^ — * vri^, 

5"*" -^ irK form the sides and B^ -^ ttt] (or 77') the diagonal of a quadrangle 

shown in FigJ^. Thus, the measurement of these 5 decay rates fixes this 

quadrangle of which the other diagonal is A3/2 (see Eqns.(^) and Fig.g) [Q]. 

^ There is a discrete ambiguity in this construction since A's ABC and ADC can be 
on the same side of the common base AC in Fig.g. 

9 
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Figure 2: The polygon construction for Method 2. Knowing IA3/2I from this 



Figure, we can determine 7 using Eqn. ([TT|) if the decay rate B^ — >• 7r°7r+ is 

also measured. 

10 



Knowing the magnitude of A3/2 from Fig.^ and the decay rate B^ -^ 11^71^, 
we can determine 7 from Eqn.(|ll]). Thus, 7 can be determined (up to a two- 
fold discrete ambiguity) by the measurement of the decay rates for 6 modes 
- 5+ — > vr'^Tr^, ttK, Bd —>■ ttK, and Bs —>■ tt^t] (or Tc^rj'). In the method of 
[0], measurement of the rates for the CP-conjugate processes of all the above 
modes is also required to determine 7. 

4 Discussions 

We comment on the accessibility of the various decay modes used in the two 
methods. The Bd and B~^ decay modes should be accessible at the e^e~ 
machines whereas the Bg -^ ttt] [t]') decay mode will only be accessible at 
a hadron machine. Since the QCD penguin does not contribute to this Bg 
decay, the decay rate is expected to be small. The measurements of the 
decay rates to CP eigenstate final states: Bg -^ nr] (r/'), Bd -^ 7r°K° and 
Bd —>■ K^rj [r]') require external tagging. 

As mentioned earlier, we have used flavor SU{3) symmetry in both the 
methods. In the factorization approximation, SU{3) breaking in the tree 
level amplitudes can be incorporated by factors of fx/ fn (see, for example, 
I]). For example, Cf^ {AS = 1) = Jk/U x Cf^ i^S = 0). However, 
since some of the strong penguin operators are {V — A) x (V + A), in the 
penguin amplitudes, the SU{3) breaking effects are difficult to estimate, but 
the breaking will still be less than ~ 0(30%). In method 2, we use the decay 
mode Bg —>■ vrr/ (r/') which does not have the QCD penguin contribution, 
but does have the EWP contribution. The EWP operators Oj^s have very 
small WCs whereas the EWP operators with significant WC's, Og^io, are 
Fierz-equivalent to the tree level operators Oi^2 [§!• So in the factorization 
approximation, the corrections due to SU{3) breaking in relating the penguin 
amplitudes for Bg — > ttt] (77') to the ones for Bd, B^ -^ ttK are given by 
factors of frj,n'/fK- 
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We have also assumed that the SU{3) breaking in the strong phases 
is small. A possible justification is that at the energies of the final state 
particles ~ 'mh/2, the phase shifts are not expected to be sensitive to the 
SU{3) breaking given by, say, mx — vn.,^ (which is much smaller than the 
final state momenta). However, it is hard to quantify this effect. 

Both the methods can be used with CP-conjugates of all the decay modes 
as well. We have neglected annihilation amplitudes: ^15, D^ and D\^. The 
validity of this assumption can be checked by comparing the decay rates 
-Bs -^ TTf] and Bg — > vrr^' - these two decay amplitudes differ only in the 
annihilation contribution (see Eqn.(^)). In the absence of significant anni- 
hilation contribution, the decay rate for Bs —>■ nr]' should be twice that for 

Bs -^ TIT]. 

In summary, we have discussed two new methods (based on fiavor SU{3) 
symmetry) to determine the weak phase 7 using the decays B^, B^ —>■ Krj 
[r]') and B^ -^ rrr] [rj'), respectively. These methods partly take into account 
rescattering effects. 
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